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Abstract. 

The operator equations for quantum hydrodynamics are discussed and solved in a 
simple cylindrical geometry. We find a solution with the velocity curl "frozen" into a 
density of the liquid in the absence of singular vortex lines. The spectrum of small 
oscillations around this solution is found as a generalization of the standard phonon 
spectrum, and the stability of the solution is demonstrated. 



1. Introduction 

With the development of physics of cooled atoms in traps, the studies of quantum 
rotation in a many-body system have revived and entered a new stage. In superfiuid 
helium-4, as it is well known [H [2], the vortices with quantized circulation emerge 
at a certain critical cranking velocity. At sufficiently fast rotation, the vortex lines 
form a lattice that on average imitates the rigid-body rotation [3]. Analogous vortices 
appear in penetration of the magnetic field into superconductors of second kind [1]. 
Roughly speaking, the first vortex can be born by a perturbation violating the angular 
momentum conservation at the angular velocity corresponding to the crossing of single- 
particle levels corresponding to the zero and non-zero angular momentum per particle. 
A nuclear analog is the so-called back-bending that occurs also at a band crossing in 
the rotating frame |5j. Similar phenomena are known for molecules [6j, predicted for 
neutron stars [7], and important for some quantum field-theoretical models [8]. 

Our study is also partly motivated by the recent explosion of activity in 
experimental and theoretical studies of superfiuidity and other quantum phenomena 
in trapped ultracold atomic systems [9l [TOl HI] and quantum Hall superfiuids [121 fT3] - 
Also, there have been many experimental studies of the fiow properties of superfiuids 
confined in porous media such as Vycor glass and containers packed with fine powder 

[a [151 [mini. 
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The question of a possibility for a quantum liquid to be in a state of rotational 
excitation without forming individual vortices or a vortex lattice is still debatable. 
Such rotation should have a moment of inertia close to that of a rigid body. A 
small quantum system, such as a medium or heavy atomic nucleus, clearly displays 
numerous rotational bands built on various intrinsic configurations. A Fermi system 
without pairing correlations should have, with semiclassical accuracy, the moment of 
inertia corresponding to the rigid-body rotation for a given density distribution [T8] . 
In agreement with the substantial presence of pairing, the nuclear moment of inertia 
is considerably smaller than that of a rigid body but still much greater than that of 
a perfect liquid [20] . Rotation of an atomic system with and without a quantized 
vortex was considered in [21] (see also references therein). The evidence for the so-called 
Landau state with identically vanishing curl of the superfluid velocity was presented in 
[22] . Especially rich spectrum of rotational phases is predicted and partially observed 
in supefiuid helium-3 [231 121], where the orbital momentum is strongly coupled with 
spin structures. New types of rotation are expected at high pressure [251 ES] and for 
crystallizing quantum liquids [27]. It is also worthwhile to mention here the fractional 
quantum Hall effect where the vorticity of Laughlin's state is proportional to the local 
density of electronic fluid [281 129] . 

The statement that a quantum liquid rotates necessarily with emergence of 
quantized vortex lines seems to be too restricting. 

In the present article we consider the quantum hydrodynamics in the original form 
developed by Landau [30] . With the exact solution of quantum operator equations of 
motion we come to the type of rotation that is not related to the formation of singular 
vortex lines (the vortex solution is possible as well). We show also that the new solution 
is stable with respect to small oscillations and find explicitly the spectrum of such 
oscillatory modes. 

2. Quantum hydrodynamics in operator formulation 

The equations of quantum hydrodynamics for a system of identical Bose-particles 
of mass m were formulated by Landau in terms of the field operators of density, 
p(x) = mn(x), and current density j(x). 



where the subscript a labels single-particle variables, and [... , ...]+ means the 
anticommutator. The quantization is performed through the commutation relations. 




(1) 



a 




(2) 




(3) 




(4) 
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and 

bi(x),Jj(x')] =i^€ijk{cnilj{x))f^5{x-x'). (5) 

The microscopic velocity field v(x) can be defined (at the points where the density 
n differs from zero) as 

, . 1 r 1 



V X = - 

^ ' 2 



-n[x.) 

It satisfies the commutation relations 



(6) 



and 



[v(x),n(x')] = -^^V.5(x-x'), (7) 



[fi(x),fj(x')] =i— eijfc^^(curlv(x))fe(5(x-x'). (8) 
m nix.) 



H 

2 



As a consequence of eq. ([7]), the density operator commutes with the curl of velocity, 

[n(x),curlv(x')] = 0. (9) 

The whole algebra can be simply postulated or derived in terms of the underlying particle 
variables, either in configuration space [301 (^^^ right hand side of eq. ([8]) is given there 
with an opposite sign), or in the secondary quantized form [3T] . 

There are statements in the literature, for example [321 ES], that quantum 
hydrodynamics formulated in such a way contains intrinsic contradictions. Indeed, with 
formal manipulations by commutation relations, one can construct the states, where the 
expectation value of n{x) is negative and very large by absolute magnitude. However, 
this is not a physical difficulty. Thinking of a real quantum liquid, we should have in 
mind the Hamiltonian of the type 

771 C 

/ d^x v(x)n(x)v(x) + W{n}, (10) 

where the functional W depends on density and its spatial gradients. Since the liquid 
consists of quantum particles, this functional contains the well known quantum potential 

W = W, + W, W,= —[d^x^^^. (11) 

8m J n 

Similarly to the centrifugal barrier in the radial Schrodinger equation, the term Wq 
does not allow the system to penetrate into the regions of n < 0. Such "fall onto the 
center" is possible only if the interaction W has a stronger singularity at n — )■ 0. We 
can conclude that the growth of quantum fluctuations at n — )■ prevents the transition 
into unphysical states. 

Being interested in the spectrum of low-lying elementary excitations in superfluid 
^He (phonons and rotons). Landau considered in Ref. [30] only irrotational motion, 
still mentioning the possibility of vortex motion in some quantum states which are to 
be separated by an energy gap from the ground state. Geilikman [31] linearized the 
equations of motion for weak excitations above the uniform background (v° = 0, n° = 
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const), found the well known Bogoliubov phonon spectrum [33], and showed that the 
modes with curl v 7^ appear only in higher orders as a result of phonon anharmonicity. 
Below we construct a class of solutions, where both, classical (c-number) and quantum 
(operator), parts of curlv are present. 

3. "Frozen" rotation 

Let us define an axial vector operator r(x) that, according to eq. (Q, commutes with 
n(x), 

curlv(x) = n(x)f(x). (12) 

Using our operator algebra, we can derive the commutation relations of the components 
Fj among themselves and with components of velocity, vj. The simplest non-trivial 
solutions compatible with the operator algebra correspond to the states of the "frozen" 
rotation [35], when F is a constant c-number vector, let say along the z-axis, while the 
density and velocity fields do not depend on z. In this case 

[curlv(x),i7] = -^fdiv[v(x)n(x)+n(x)v(x)] =z^f^^. (13) 

Therefore, for time-independent density, the quantity curl v(x) is the constant of 
motion that characterizes the stationary states of the system. 

We consider a planar motion v(r) in a circular cylinder of radius R with longitudinal 
axis z, Vz = 0. The transverse components of the velocity field satisfy, according to eqs. 
(IH]) and (fT2|) . the commutation relations similar to those in a magnetic field, 

Kir),vy{r')]= ^r5(r-r'), (14) 

where r = (r, (p) is the coordinate vector in the (x?/)-plane. The operator equations of 
motion for the velocity field, on this class of states, are given by 

r / \ rn ih d fmf^(r) 6W ] ..^ 1 r / x / M /-,^\ 

Here the density n(r) has to be expressed in terms of curlv due to eq. (!T2l) . 

We will look for the solution of the set f fT^][T5l) as small oscillations above 
the background characterized by the classical fields of velocity, v°(r), and density, 
determined by curl^ v°(r) = Fn°(r). We limit ourselves by the axially symmetric case, 

v;{r) = u{r), < = 0, n°{r) = u{r). (16) 

Of course, the angular momentum of the system will be directed along the axis but, due 
to the azimuthal symmetry, there is no need to make a transformation to a rotating frame 
and introduce the Lagrangian multiplier of the angular velocity. The similar situation 
is known also in nuclear physics where the rotational motion around the symmetry axis 
does not require the cranking procedure. The difference here is that we find collective 
excited states instead of the more or less random sequence of single-particle excitations 
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carrying the angular momentum projection onto the symmetry axis in a nucleus rotating 
around the same axis. 

The background state, as follows from eq. ( 1T5|) . should satisfy the equilibrium 
condition, 

mw^ _ d f6W^ 
r dr 



5n 



[17) 



n=v{r) 



In our geometry. 



du u ^ , . 
— = - - + Vv(r). 
dr r 

This leads to the general form of the background velocity field. 



u{r) 



K + T I dr' r'v{r') 





K+ -N{r) 
2,11 



(19) 



where is a constant and N{r) is a number of particles inside the radius r, < r < i?; 
here and later all macroscopic quantities are taken per unit length of a container (see 
also Appendix). 

Our general solution (1191) represents a superposition of a vortex line of intensity 
K and rigid-body rotation of cylindrical layers with a local angular velocity VL[r) = 
ri/(r)/2. The angular momentum of the liquid L = can be found as 

L = m J d^rv{r)ru{r). (20) 

Using eq. (fT9l) . we derive 

"^iv^ (21) 



L = mNK + 



where = N{R) is the total particle number. In the absence of F, as usual, the bosonic 
quantization L = Nhs, where s is an integer, determines Ks = sh/m. The additional 
part of the angular momentum does not carry any singularity at r — )■ 0. 
For a vortex- free case, if = 0, we have from L = Nhs 
Anh 



mN 



(22) 



which corresponds to rotation with local linear velocity Vb{r) = Fz/(r)r/2. 
4. Solutions with and without a vortex 



To move further, we assume for simplicity that the functional W does not contain 
gradients of n and can be expanded around the unperturbed density of non-rotating 
liquid, no. 



W = W{no} + ^m'^g / d'^r[n{r) - Hq]'^ 



(23) 



With this functional, the balance equation f lT7|) takes the form (see Appendix for 
derivation) 



mu 



2 

m g— + 
dr 



h^ d 


1 


(rfr) 


2 d 


8m dr 


z/2 




ur dr 



dv^ 
dr 



(24) 
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The compressibility g defines the speed of sound, c = ^Jfngn^. It is convenient to 
introduce dimensionless variables and express velocities in units of c, energies in units 
mc^, lengths in units of tq = h/mc, and density in units of no- Then the only remaining 
parameter is V [in usual units V =^ huoT/ (mc^)]. For typical quantum liquids, the 
length ro corresponds to atomic distances, while the macroscopic container has 1. 
At the same time, for reasonable angular velocity, the rotation is subsonic, TR < 1, 
which means F <^ 1 in reduced units. 

As follows from eq. ( lT9l) . we can have two types of solutions for eq. ( l24l) with a 
regular behavior of density i/(r) at the axis, r — )■ 0. For the vortex line solution, K ^ 0, 
at small distances from the axis the velocity is u ~ K/r, while the density goes to zero, 
z/(r) oc r^l^l + (9(r^''^l"'"^), and its behavior is determined by the quantum term ( ITTj) . 
Here F influences only higher terms of the expansion of i^(r), the radius of the vortex 
core does not noticeably depend on F and has the magnitude of the order of one (the 
distance where we can equate contributions to the energy from compressibility and from 
quantum effects). Such results for quantized vortices are well known [H |36]; they are 
modified by rotation only at large distances from the center. Below we will be interested 
in the solutions of the second type, where the vortex singularity is absent, K = 0. 

At K = 0, the velocity and density satisfying eqs. and (fT9|) can be expressed 
by a well converging series in powers of a parameter ^ = uqT'^R^: 

Hr)-o|l+|^ + ^^ + ...}. (25) 

u(r) = - z/oFr ll + —— + j + ...} ■ (26) 

^ ' 2 ° \ 16 R^ 384 R^ j ^ ^ 

The central density z/(0) = is determined by the particle number conservation being 
related to the unperturbed value uq = 2u{R) / {RT) = 1 through a similar series, 

= 1+ — + — + ... (27) 
° 16 384 ^ ' 

It is easy to estimate that in such a state the contributions to energy per particle 

from the kinetic term and from compressibility are, by order of magnitude, ~ T^R^ and 

~ F^i?^, respectively, whereas the quantum potential brings in only a small contribution 

~ F^. The angular momentum, eq. f l2ip . equals L = Nhs, where s = TR'^/A; at TR ~ 1, 

this value s ^ 1. In the case of the external cranking with angular velocity fl, the state 

with a certain value F 7^ is getting energetically favorable at a critical value 

For the minimum value Fmin = F(s = 1) we have ^ R^^. This type of 

rotation can emerge at a smaller (by a factor of In R) angular velocity than the creation 
of usual quantized vortices. However, the probability of such an excitation, for example 
by an accidental nonuniformity of the walls of the container is considerably lower, 
while the observation at TR < 1 is difficult. 
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5. Small amplitude oscillations 

Now we will construct the low-lying excited states corresponding to small oscillations 
of the density and transverse velocity around the classical solution [eqs. fl25ll2l7)]. We 
introduce the exciton operators, v(r) = v — u and n(r) = n — v. The linearization of 
the equations of motion ( |T5i) in the model ( l23i) gives 

[v, H] = -zV(u ■ V + P) + i[{uh + iru) x f], (29) 

or, in cylindrical coordinates, 

. d 



and 



-I Tri^^v + P) + ^r({)^z/ + un) 
or 



r OLp 



iVvj-v. 



Here the operator was introduced. 



P 



n 



n 



t VM- -— -V- 




and, according to eq. (IT2 



Tn 



1 

r 



d 



dVr 



dr dip 



The continuity equation. 



[n,H] 



u dn 



d 



dr 



^ruVr) + iy 



dip 



1 

r dip r 

follows from the equations of motion. Finally, with use of eq. (H 
the equations of motion in a symmetric form, 

. dP 



[Vr,H\ + I - — 2l-V^ 

r dip r 



dr 



u dv^ ^ „ i dP 

r dip r dip 



(30) 
(31) 

(32) 

(33) 
(34) 

we can write down 

(35) 
(36) 



At azimuthal symmetry of the ground state, the elementary Bose-excitations can 
be characterized by a certain angular momentum = M. The quantization goes with 
expansion of our fields over corresponding creation, AIi, and annihilation, A^i, operators, 
where the additional subscript k labels consecutive modes for given i. The expansion 
can be written as 



V, 



.(r \ 



/ 



n(r) 



E 



vlf-^r) 



n 



E 

ke 



{RM{r)AM + Rl^e{r)AU} 
i[Fk,{r)Ak,-FU{r)Al^,} 
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Here the unknown radial form-factors of excitations, Rke, ^ke, and F^i, are the matrix 
elements between the ground state and excited one-quantum state with quantum 
numbers k, £ and energy Eke- They satisfy the set of linear differential equations, 
„ „ M \ „ u dPuf 

Eki-i-)Rki + 2- ^ki = -7^, 38 
r / r dr 

Eke-i-) $H + TuRki = - Pke, (39) 
r / r 

m, = — ^ + — - -Pfc^ (40) 
dr r r 



The right hand side quantities Pkt are related to the matrix elements of density: 

Fkt 



Pki — Fke — - 



^2 1 dv d 

^ V dr dr 



5 



(41) 



where the differential operator V| is defined as 

_2 Id 

dr'^ r dr r 

To analyze this set of equations, it is convenient to introduce the notations 

£ki{r) = Eke-i -, Dkiir) = £l - 2T — , (43) 

where Ske{r) can be interpreted as "excitation energy" in the frame rotating together 
with the layer r = const. The matrix elements of the velocity components can be 
presented as 

Rk,{r) = D^l [Sk^PU - 2£ ^ P,,) , (44) 

and 

^ki{r) = D,l {Sk, ^ Pki - TvP',)i ■ (45) 
here and below the prime indicates the radial derivative d/dr. Returning to the definition 



0]), we obtain a differential equation for Fkt. 



DktFkt + u'V'jPki 



„2 



As seen from eq. ( BTj) . the right hand side of ( H6|l is the third order differential operator 

acting on Fki] the coefficients of the operator contain only small terms of the expansions 

( l25l) and ( l26ll . Therefore we can act iteratively. 

Here we will limit ourselves by the main order assuming 

u T 1 
u ^ uo^l, ~ ^ 2' ^ki ^ Eke - -£T = const. (47) 

Then eq. fH6|) is reduced to 

{4e-r' + V'e-\^t)Flf{r) = 0. (48) 
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The regular at the origin solution is a superposition of cylindrical functions 

Fi?{r) = ht^Uq^^\) + (49) 
where the wave numbers g^^^ are related to energy E^e, by a dispersion relation 

or, selecting the appropriate roots, 



i?' = 2(vr+^i^Ti). (51) 

In the same approximation we obtain for the function fj4T|) : 



\ {t^iTjM^M^r) - fet^git^'/KgtV)} . (52) 

We came to the direct generalization of the standard phonon spectrum for the case of 
the frozen vortex-like motion. There is now a small gap F in the phonon spectrum of 
transverse oscillations. The gap is determined by the magnitude of curl v. 

The coefficients fc^^'' of the superposition ( l49l) . as well as quantization of the 
eigenvalues Em-, are determined by the normalization of the solutions and the boundary 
conditions at r = i?. For our problem of vibrations on the background of a classical 
solution with almost constant density, the natural boundary conditions in the accepted 
approximation are those of vanishing normal components of the velocity, Rki-, and of the 
gradient of density, F^^. A detailed discussion of boundary conditions can be found in 
Ref. [38]. With the aid of eqs. (I41|49l) and fl52l) . we find, up to a normalization constant 
Cm, 

= CMl'AqM^R), bi-^ = -CMJ'eiqlVR)- (53) 
The energy quantization is given by {z± = q^^ R) 

As it follows from eq. (I5T1) that g*^"-*^ — g'-^^^ = 4, we need g*-^-* > 1. For a macroscopic 
vessel, R^ 1, and i < R, the cylindrical functions li and in the last equation can be 
taken in asymptotics. This leads to 

^Mqte\ql:^" + q^^'m^^) - r I {qt^^'Uz^) + git^') = 0. (55) 

In the irrotational case, F = 0, or for £ = 0, the dispersion equation flS^ determines 
the usual hydrodynamic spectrum z = zl^, where J'i{zlg) = 0. Assuming a small shift 
6zke of eigenvalues, we find (omitting subscripts k, i) 

5z ^ AR' + jzr I ^ ~ i 

z° ^ 4R^ + 2{z°y\ ^ + {z°y + {z°)y4R^ p + {z°y' ^ ^ 

where the parameter ^ was used earlier, eq. fl25|) . This result is valid for short 
wavelengths, when g^"*") ~ g'^^^ ~ 1, 2;° ~ i? ~ and 5z/z° ~ T/i <^ 1, but also 
in the long wavelength limit, g^^-* ^ 1, if £ is sufficiently large, 1 <^ £ ~ q^^^'R <^ R. In 
the last case it is applicable even for ^ ~ (^°)^ since 5z/z° ~ <^ 1. For excitations 
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with small values of quantum numbers ~ 1 and F ~ 1/R, the displacements of the 
roots is significant, and one has to return to the general eq. fl5^ . In that case, T ~ q^^\ 
and the presence of the gap in the spectrum ( 150|) is essential. 

The simplest way to finding the normalization constants is based on the 
orthogonality conditions following from the set of equations ( l38lH0l) . 



R 



-Lcfi)Wcfi)(°) ip{0)p(0) \_ 1 =fc 1 Q r 



kiOkk' 



(57) 



where the right hand side quantities Bfc^, after some algebra using the equations and 
boundary conditions, can be expressed as 



e 



Ti 



ke 



DuS-ki 



+ 



D 



ke JO 



R 



drrFtVPi 



(0) p(0) 
ke ^kl ■ 



On the other hand, conditions fl57j) and operator expansions fl371) show that 



A 



ke 



Qke Jo 



(5^ 



(59) 



The quantization according to the Bose-statistics, [Ake,Al,^,] = 6kk'See', together 
with the original rules of quantum hydrodynamics ([7]) transformed in our problem to 
eq. f|T^ . leads in this approximation to 



92 



ke 



[drr[VRtM} + £4(R'SY + (4?)^ } = 1 



The explicit calculation of the integral fpU]) and comparison to eq. 
simple normalization condition, 

Qke = 2£^fc£. 



(60) 

determines the 
(61) 



The same result can be derived by the direct diagonalization of the linearized 
Hamiltonian. We will not write down a cumbersome expression for the normalization 
constants Cke that can be derived from eqs. (!58|) and f l6T|) with the aid of the solutions 
[9|) and ( 152|) . Here we show only the result for a pure hydrodynamic spectrum that 



becomes exact for 



0, when b^^^ — )■ 0, while 



(+)^2 
ke ) 



8D 



ke 



(qie^fJeiQlVR) 



SkfR 



(?l?^)V 

In spite of the presence of the gap in the energy spectrum (l50ll . we still obtain in 
the transition to an infinite system, i? — > oo, a vanishing interaction strength with the 
long-wavelength, q^^^ — )■ 0, phonons. Indeed, in agreement with general arguments, we 
see from eq. (162|) that 

^ ' ^ (63) 



.(62) 



^ke 



Dke oc q^^\ 



It is easy to see also that the off-diagonal matrix elements of the local density are small 
being at F ~ q^'^^ ~ 1/i? of the order N~^^'^ as it should be for collective excitations. 



(H|n(r)|0) 



kl ^' 



Je 



lie 



(64) 
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This justifies the hnearization similar to the random phase approximation. 

Our energy spectrum (150|) formally allows for two types of excitations with phonon 
angular momentum i and energies 



In the long wavelength limit, the Coriolis interaction, —£^l, of the "minus" phonon with 
^ > 2 is so strong that the creation of such a phonon becomes energetically favorable 
at sufficiently large angular velocity, namely if F > 2{[(g^,^^)^ + (q'it^)''/4]/(£^ — 4)}. 
However, this would correspond to the supersonic velocity of rotation at the boundary, 
TR/2 > 1, when the approximation used above is invalid and the spatial non- uniformity 
is essential. In our approximation we have to take into account only the "plus" branch 



6. Discussion 

The study of rotating helium II has a long history. Contrary to the predictions of Landau 
theory [20], it appeared that below the A-point the superfluid components behaved in 
rotation as a classical fluid. This was shown by measuring the shape of the free surface 
meniscus of rotating helium that revealed a classical profile [391 HHl SH HSl SSj HI]. 
The measurement of the angular momentum of rotating helium II was also found to 
correspond to the rigid body rotation of the liquid[45l HHl H?]. As is well known, this 
paradox was solved by the Onsager-Feynman theory of quantized vortices [3]. The critical 
velocity for the vortex nucleation is Qg = {fi/mR'^)ln{R/a), where R is the radius of 
the vessel and a is the radius of the vortex core. In the aforementioned experiments the 
speed of rotation was much higher than Qg so that the vortex lines formed a lattice that 
on average imitated the rigid-body rotation [3]. 

On the other hand, we cannot exclude the manifestation of "rigid-body" rotation 
in several reported experiments. For example, the interpretation of some data in [22j as 
Landau state does not seem unambiguous. In this paper the helium was contained in a 
torus resonator filled with slip rings and brushes so that the normal-fiuid component in 
rotation was locked to the resonator by its viscous interaction. The relative velocity of 
the normal fiuid and superfiuid was determined by measuring the frequency difference of 
fourth-sound modes which run around the resonator in both directions. The dependence 
of the superfluid angular velocity on the normal-fluid angular velocity had the form of 
a hysteretic curve. The reversible part of this curve was attributed to the vortex-free 
potential flow (Landau state). Plausibly, this reversible part could be as well attributed 
to the "rigid" velocity fleld, since for the multiple-connected geometry the contribution 
from such a rotation of superfluid (if it does exist) cannot be distinguished from the 
contribution of irrotational potential motion. 

It is also appropriate to mention here the paper [18], where the angular momentum 
of superfluid helium was measured after the rotating helium bucket has been stopped. 
The resulting angular momentum was smaller than the angular momentum of rotating 




(65) 
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bucket, and the authors came to a conclusion that the reason for this was the formation of 
the long-hved macroscopic rotational state that accepted part of the angular momentum. 
As the authors pointed out, this fact could not be explained neither by the Landau 
two-fluid theory (in a simply connected region, the only solution for potential flow 
of incompressible fluid is f = 0) nor by Onsager-Feynman vortex states. In this 
connection it is worthwhile to remember the so called Lin's hydro dynamics |^ . where, 
to our knowledge, for the first time the conception of rotational superfiuid flow has been 
proposed. Lin has suggested that the superfiuid might have very small nonzero viscosity 
and would thus reach a steady state of uniform rotation in which boundary slip would 
cause the superfiuid to rotate more slowly than the normal fluid at low velocities. Lin 
has applied his approach for the explanation of the variety of helium-II experiments. In 
particular, the above mentioned experiment [48] has found a satisfactory explanation 
within his theory [50]. 

For our velocity field ( (fT9|) with K = 0), it is easy to calculate critical velocities for 
slow rotation assuming that the density is constant over the volume. For the rotation in a 
circular cylinder of radius R we obtain, by minimizing the free energy in a rotating frame, 
the critical angular velocities for vortex-free rotation Qi, = ^o(2p+l), where p = 0,1, 2.... 
Here ^Iq = h/mR^, and (see ([22])) s = ioi Vt < VLq, s = I ioi VLq < Vt < 3i7o, s = 2 for 
3f2o <VL < 5i7o 3^iid so on. Comparing this with the critical velocity for the appearance 
of the first vortex Vlg, we see that the first critical velocity for vortex-free rotation VLq 
is approximately ten times smaller than the vortex critical velocity VLg. In general, our 
critical velocities are on the order of VLq which even for R \ mm are rather small, 
about 0.01 rad/s. 

Now that the instrumentation methods are greatly advanced as compared to those 
used in early experiments, a possible clear-cut experiment for a search of the frozen 
rotation would be the measuring of the free surface of rotating helium with the speed 
of rotation being much less than the critical velocity for the vortex nucleation. Other 
possibility is the direct investigation of the flow pattern below the A-point with the aid 
of some visualization technique (see[5T] and references therein). It is also worthwhile 
to note the paper [52] where the circularly polarized third-sound wave modes were used 
to drive persistent flow states in a superfiuid helium film inside the resonator. As was 
pointed in this paper, this method could, in principle, be used to find the profile of the 
velocity field. 

One can also point out a possible connection of our solution to the effects emerging 
in the flow of helium-II in porous glass with the size 1-10 nm of the pores [53], narrow 
slits [54], in small tubes [55], and submicron channels [56] |57]. In particular, as was 
reported in [56], the critical velocities for the flow of superfiuid ^He through a submicron 
aperture appeared to be much smaller than it would be necessary for a vortex ring 
nucleated at the wall. 

An interesting question of experimental manifestations of a "rigid-body" rotation 
of quantum Bose-fiuids other than liquid helium requires special consideration and is 
beyond the scope of the present paper. 
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In this article we made an attempt to vindicate the old operator formulation of quantum 
hydrodynamics due to Landau [3_0j and Geilikman [31]. This formulation does not lead 
to any unphysical states of the type of the "fall onto center" and allows one, in the 
simplest cylindrical geometry, to write down and solve the operator equations of motion 
for a quantum liquid. An example of the explicit solution was found in the form of the 
classical background and collective oscillatory modes. The solution in general contains 
a superposition of a vortex line of intensity K that should be quantized through its 
relation to the angular momentum and a rigid-body rotation with a local angular velocity 
Q{r) = ri/(r)/2 due to the velocity curl "frozen" into the density of the liquid. 

We considered in detail the solutions without a vortex singularity, K = 0. 
The calculations are rather long but straightforward. On top of such a solution we 
superimpose small oscillations of transverse components of the velocity and related 
oscillations of density and linearize the equations of motion. For the spatial function 
we get a complete set of modes with orbital momentum i and the radial dependence 
determined by the combination of cylindrical functions Je{kr) and Ie{kr) that is 
determined by the boundary conditions. 

Now, when the new types of quantum liquids are under experimental scrutiny, the 
search for diverse manifestations of rotational superfiuid flow is at the frontiers of cur- 
rent interest. In any case, it seems that the conventional consideration of the vortices 
as the only long lived objects in rotating quantum liquids is too restrictive. 

We acknowledge useful discussions with S.T. Belyaev. We also thank the referee 
who pointed out the similarity between our frozen rotation and the vorticity of 
Laughlin's state in fractional quantum Hall effect. 

The research was partly supported by the NSF grant PHY-0758099. 
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To find the ground state based on the rotational velocity profile, 

u{r)=^N{r), (66) 
/vrr 

we have to minimize the functional depending on z/(r) and its gradient dv/dr. The 
functional contains kinetic energy /C, potential (interaction) energy W, and quantum 
potential Q. Let us consider these terms one by one. 
The kinetic term is 



JC = — J d'rpu\ (67) 

Here we minimize directly taking into account that we restrict ourselves by the class of 
states, where eq. ( 1^ is valid, 

S}C= \5iy(r)N\r) + 2i^(r)N(r)6N(r)} . (68) 

Here, using the step function 6(x), 

m{r) = [ d^r'e{r - r')5v{r'), (69) 



and changing notation in the double integral r ^ r'. 

The potential term, 

W=^m'gJ dh{v-n,)\ (71) 

is trivial, 

d'^r{u -no)6u{r). (72) 

The quantum term, 

Q^fl^r^^. (73, 

8m J V 

depends both on v and Vz/. The variation goes as in the Lagrangian mechanics, 

^"2- |l/dV{§^..+ ^.,(V.)).(74) 



Here 



dv 
dC 2 



(75) 



(Vz/-5(Vi/)). (76) 



diyu) V 

Since the symbols b and V commute, eq. f l7^ gives 

/Vz/ \ / Vz/\ /Vz/ \ V^z/ (Vz/)^ , , 

2V- ( — 5z/) -2(V- — ) =2V- ( — 5z/) -2-— +2^-^. (77) 



Z//\Z// \ V J V V 
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Combining eqs. and (17^ . we obtain 



where we have also a surface term, 

If we forget about the surface contribution, the Lagrange equation is 

m^^ zz-no + — V'—f- -2 80 

8m [ z/^ V j 

mV^ {N^ir) f<Pr' ^ 



+ 2 / ^ 9{r' - r)u{r')N{r') \ . 



Taking the gradient d/dr in eq. (15U]1 and using again eq. (jHS]), we come to eq. 



